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Abstract 

Quantum Monte Carlo method is applied to fractional quantum Hall systems. The use of the linear programming method 
enables us to avoid the negative-sign problem in the Quantum Monte Carlo calculations. The formulation of this method 
and the technique for avoiding the sign problem are described. Some numerical results on static physical quantities are also 
reported. 
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1. Introduction 

Two-dimensional electron gas (2DEG) in a strong 
magnetic field is known as a system wliicli sliows vari- 
ous plienomena resulting from electronic interactions. 
In particular, the fractional quantum Hall (FQH) ef- 
fect has attracted much tlieoretical and experimental 
attention as a correlation-dominated one peculiar to 
tliis system. To investigate such strongly correlated 
systems theoretically, numerical studies play a signifi- 
cant role. For example, the exact diagonalization and 
density matrix renormalization group (DMRG) meth- 
ods [1] have been used for the study of FQH systems, 
because the ones are directly accessible to the ground 
state. However, the exact diagonalization can study 
static and dynamical properties only for small systems, 
while it is not so easy to obtain dynamical information 



by the DMRG calculations. Thus we apply the quan- 
tum Monte Carlo (QMC) method to this system in the 
present work. 

The QMC method has been used for several quan- 
tum many-body systems. It can investigate both static 
and dynamical properties of systems larger than those 
to which the exact diagonalization is applicable. How- 
ever, this method is inherently involved with the so- 
called negative-sign problem. Thus we need to over- 
come (at least, or moderate) the negative-sign nuisance 
in order to make the best use of the method. In this 
study, we formulate the QMC method to the FQH sys- 
tem that is free of the sign problem and report some 
numerical results. 
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2. General formalism of QMC at zero 
temperature 

In the zero-temperature formalism, the expectation 
value of an observable O for a Hamiltonian H is given 

by 

{0)= Mm. {'4)\e-'^^^^Oe-'^^'^\'4>)l{'i)\e-^^\'4>), (1) 

where \'tp) is an arbitrary state not orthogonal to 
the ground state. Here we assume that the Hamil- 
tonian can be written in a quadratic form of a one- 
body operator. When the imaginary-time evolution 



operator e 



-I3T-L 



is decomposed to imaginary-time 



slices e~^^^ • ■ • e~^^^ , we perform the Hubbard- 
Stratonovich (HS) transformation for the quadratic 
form by introducing auxiliary fields for each time slice. 
Then the Hamiltonian is reduced to a linear form of 
the one-body operator and the expectation value is 
expressed in a form of an auxiliary-field path integral. 
We evaluate it by means of the Monte Carlo method. 
However, statistical weights used in this Monte Carlo 
calculation are not always positive for any auxiliary- 
field configuration and then the normalization denom- 
inator in Eq.(l) often becomes almost zero. Thus it is 
often impossible to calculate (C) with sufficient preci- 
sion. This disaster is called the negative-sign problem. 



3. Hamiltonian 

We consider interacting electrons confined on a 
spherical surface with a magnetic monopole located at 
the center of the sphere [2] . For simplicity, we neglect 
the spin degrees of freedom and assume that all the 
electrons occupy the lowest Landau level. Then single- 
particle states are specified by the z-component, m, of 
angular momentum whose amplitude is s, where 2s is 
the number of flux quanta piercing the sphere and m 
ranges from —s to s. The Hamiltonian can be written 
in a quadratic form of the density operator and its 
time reversal: 



^ 2s K 

■9 PknPkn +GoPaa, (2) 



Pkn = E {KN\rmm2)i-l) 



where am is an annihilation operator of electron, 

{KN\inim2} the Clebsch-Gordan coefficient, and Co is 
a constant. The time reversal of the density operator is 
defined by p^^.^ = {—Vj^p^,^. The coupling constant 
Xk is related with the Haldane pseudopotential Vj as 

2s 

Xif = ^ TkjVj, (3) 



= (-1) 



{2J + 1) 



s s J 



s s K 



where the braces denote Wigner's 6j symbol. Since the 
transformation matrix T satisfies = T, we also 
have a relation of 

2s 

Vj=Y^ TjKXK. (4) 

x=o 

The Hamiltonian in Eq.(2) can be written in a linear 
form of the density operator by the HS transformation. 

Denoting the auxiliary field for a mode (KN) by (Tkn, 
the linearized Hamiltonian is given by 



(<y*KNPKN + CTknPKN 
1 



+ 2X0P00 + Copoo, 



where tik = i when xk < and tik ~ 1 for non- 
negative Xk- We note that poo is proportional to the 
number operator of electrons. Since the number of elec- 
trons is a conserved quantity in the zero-temperature 
formalism, the term XoPooPoo and Co Poo in Eq.(2) can 
be considered as constants. 



4. Sign problem 

In order to discuss the sign problem, we formulate 
here the auxiliary field path integral in a matrix rep- 
resentation. We first define the matrix elements of the 
linearized Hamiltonian by 

s 

-Al3h{a) = Mm,ni<r)ai,an, (5) 



m..n= — s 



where 



am {m > 0), 
{-iy-"'am {m < 0). 
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Next, we represent the state \tp) in Eq.(l) in terms of 
a matrix V as 



where is the number of electrons. We remark that 

M is a square matrix of dimension 2s + 1 and V is a 
(2s+l) xiVe rectangular matrix. Then the denominator 
of the Eq.(l) is written as [3] 



i K=l N=0 



i Kk=1 



^P\XKo\ f A/3\xkn\ 



4tt 



X exp 



(i) |2 



= det 



(6) 



where A't = (5/ A/3 is the number of the imaginary-time 
slices. 

The sign problem is brought about by the fact that 
C(cr) can be negative. However ^(cr) is forced to be 
non- negative under the following conditions [4] : 

(i) 2s is odd and A^e is even, 

(ii) the coupling constants in the Hamiltonian satisfy 
Xif > Ofor /r = 1,2,--- ,2s. 

We show this as follows. From the condition (i), 
the indices, m and n, of the matrix M{a) in 
Eq.(5) take the values, ±^,±|,-- - ,±s. The matrix 
M{a) under the above conditions satisfies two rela- 
tions as M_,„,_„((j) = Mj^_„(cr) and M-m,nicr) = 



satisfies VL 



m,Nc/2 + i — I'm.i 



= TCi and F_m,i = -K 



m,JVe/2 + i 



for positive m and i = 1, • • • ,iVe/2, then the matrix 
in Eq. (6) can be written as 



ytgM(a("t))...gM(a(l))^^ 



A B 

-B* A* 



where A and B are square matrices of dimension iVe/2. 
This type of matrix always has paired eigenvalues that 



are complex-conjugate each other. In fact, denoting a 

matrix of the above form and its eigenvalue by X and 
A, respectively, the complex conjugation of the proper 
equation yields the one on A* : 



= det 





-1 




1 








X 




-A* 




1 




-1 





= det [X - A* 



That is. A* is also an eigenvalue of X and thus the de- 
terminant of this type of matrix is always non- negative, 
which leads to the desired condition that C,{(t) > 0. 

Now we consider which FQH systems satisfy the 
above conditions. The condition (i) can be satisfied, for 
example, in case of 2s = SA^e — 3 for the u — 1/3 Laugh- 
lin state [2], or 2s = 2Ne — 3 for the Pfaffian state [5]. 
However, the condition (ii) is not satisfied when the 
values for the Coulomb interaction are used for Vj in 
Eq.(3). Thus we need to control the value of xk (that 
is, that of Vj) by solving a linear programming prob- 
lem in order to satisfy the condition (ii) . Namely, tak- 
ing into account that only Vj for odd 2s — J are phys- 
ical for fermionic systems, we inquire which set of xk 
minimizes 

2s 

F{Xo,Xi,--- ,X2s) = X] >^j{'^TjKXK -Vj) 

2s-J:odd K=0 

(7) 

under the conditions: xo ^ Oj Xk > for K — 
1, 2, • • • , 2s, and Xj{J2k=o TjkXk -Vj)>0 for odd 
2s — J. We note here that Aj controls the variance of 
pseudopotential from that for the Coulomb interaction 
for each J. In the present study we always set \j = 1. 
Negative xo is needed to make Vb (= ^k'^okXk) 
positive under the fact that Tqk < for all K. Thus 
obtained xk minimizes the variance from that for the 
Coulomb interaction satisfying the condition (ii) . 



5. Numerical results 

Figure 1 (a) shows an example of pseudopotentials 
free of the negative-sign problem. Although the opti- 
mized potentials do not coincide completely with those 
for the Coulomb interaction, their monotonical depen- 
dence on J is realized naturally. We show another ex- 
ample of sign-problem free potential in Fig. 1(b). To 
get the short-range components closer to those for the 
Coulomb potential, we bring long-range components 
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Fig. 1. Pseudopotentials for the Coulomb interaction and 
negative-sign-free potentials obtained by the linear pro- 
gramming method for 2s = 21. The referred potential in 
the linear programming method is the Coulomb potential 
in (a), but in (b) long-range components, Vq and of 
the referred potential are set to be zero. Only the physical 
(2s — J: odd) components are shown. 
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Fig. 2. Expectation values of (a) energy per particle and 
(b) total angular momentum by the QMC calculation are 
shown versus inverse temperature for 2s = 21 and = 8. 
Optimized potential in Fig. 1(b) was used. The values used 
as imaginary-time slice width are A/3 = 1/8, 1/16, and 1/32 
in units of l/{e'^ /Ig). The ground state energy obtained by 
the exact diagonalization is shown by an arrow on the right 
side of (a). 

close to zero. Namely, we set Vo = V2 = in Eq.(7). 
As a result, variance of short-range components from 
those of the Coulomb potential is squeezed. We note 
that optimized pseudopotentials can be obtained with 
less variance for higher Landau levels or in the presence 
of finite-thickness effects. 

In Figure 2, we show expectation values of energy 
and total angular momentum against the inverse tem- 
perature /3. Although finite value of A/9 gives rise to nu- 
merical errors in the Suzuki- Trotter decomposition, ex- 
pectation values almost saturate for A/3 < 1/8 {e^/ls) 
{Ib = \fc[eB: the magnetic length). For comparison, 
we show the ground state energy obtained by the ex- 
act diagonalization by an arrow in Fig. 2(a). It is noted 
that the total angular momentum of the true ground 
state is zero. 

We find that for /3 > 20/ (e^ /Ib) these two quan- 



tities converge to the values expected for the ground 
state. In fact, the value of energy gap is about O.le^/is 
in our exact diagonalization study, and thus the (3- 
value choice of /? > 20 is appropriate to investigate the 
ground state properties. In case of smaller energy gap, 
larger values of /3 are needed for the convergence of 
expectation values. Thus incompressible gapful states 
are more suitable than gapless ones for the study by 
the present method. 

It is significant to discuss whether the ground state 
for negative-sign-free interaction is similar to that for 
the Coulomb interaction. Unfortunately the overlap 
between these two states for u = 1/3 is not so large 
(~ 0.38) for the pseudopotentials shown in Fig. 1(b). 
Since the negative-sign- free interaction with V2S-1 
quite larger than V2S-3 is expected to stabilize the 
Laughlin state as the lowest-energy state, the small 
overlap value seems to come from the reduction in 
long-range components. However, taking into account 
the fact that the ground state is gapful and the overlap 
is finite, the ground state for the negative-sign-free in- 
teraction is considered to capture essential properties 
of the Laughlin state. 
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